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^-j Abstract. We study the properties of unbound dark matter particles, defined as having 

^ sufficient kinetic energy to escape the boundary of their parent halos. We find that the mass 

i_i fraction of unbound particles increases strongly towards the edges of halos and in denser 

^ environments. Mergers are the largest source of unbound particles; however, major mergers 

are not especially more efficient than minor mergers in this regard. Because unbound dark 
matter particles can leave halos more easily than their baryonic counterparts, we also find 
that standard ACDM cosmology predicts a redshift-dependent boost for the baryon fraction 
in massive clusters which reaches as high as an extra ~ 20 percent at z = compared to the 
QQ universal value. Our findings have relevance for dark matter detection experiments, precision 

(O calibrations of the halo mass function, the use of baryon fractions to constrain dark energy, 

and searches for intergalactic supernovae. 
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1 Introduction 

Unbound particles represent a small fraction of the mass in dark matter halos [1], but they 
carry the highest energies in the dark matter velocity distribution function. This property 
means that unbound particles can have unexpected importance for experimental and ob- 
servational interpretations, including those for dark matter detection experiments, precision 
cosmology surveys, baryon fraction measurements, and searches for intergalactic supernovae. 

Direct detection experiments are sensitive to the velocity distribution of dark matter. 
This has led to many studies of the expected phase-space distribution of dark matter particles 
[2-10]. Although conclusive results have yet to emerge from the experiments themselves, 
many groups have recently started to probe the range of collisional cross-sections and particle 
masses most favored by theoretical physicists [11-15]. With possible discovery on the horizon, 
it is important to better constrain the velocity distribution of the Milky Way, as well as to 
better understand the sources of halo-halo variation in the high-energy (e.g., unbound) dark 
matter fraction. 

Another reason for studying unbound particles comes from observational surveys which 
aim to place precision constraints on cosmology (e.g., BOSS, DES, BigBOSS, Pan-STARRS, 
eRosita, Herschel, Planck, JWST, and LSST; [16-24]). In order to fully realize their statistical 
power, these surveys depend on precision calibrations of the halo mass function to 1% or 
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better accuracy [25, 26]. Traditionally, halo mass functions have been calibrated from gravity- 
only ("dark matter") simulations [27-31]. These calibrations are expected to be inaccurate 
within several galaxy scale lengths due to baryonic feedback effects [32, 33], but there is also 
another bias at large scales (e.g., at the virial radius) due to the effects of unbound dark 
matter particles. 

Specifically, high-energy dark matter can freely enter and leave halos, but high-energy 
baryons will collide with the gas already present in the halo and remain behind. Dark matter 
simulations capture the first effect but not the second, leading to an underestimate of both the 
halo density at large scales and the halo mass. This effect also suggests that the local baryon 
fraction should increase beyond the cosmic average on the outskirts of clusters. In fact, this 
has already been detected in simulations [34] and in observations [35]. However, because 
dumpiness in the intracluster medium can also produce a similar observational signature, 
proper interpretation of the result in [35] requires exploration of the relative contribution 
from both effects. 

Observations of intergalactic supernovae [36-38] also relate to unbound particles, es- 
pecially observations of supernovae traveling at high speeds [39]. Such high- velocity super- 
novae might either come from the intracluster stellar population [38, 40], from merging events 
[41, 42], or from hypervelocity stars ejected from the central galaxy [43]. Since stars are ef- 
fectively coUisionless and so have similar dynamics to dark matter, one can place an effective 
upper limit on the expected high-velocity supernova rate from stars pre-existing in the in- 
tracluster stellar population by analyzing the unbound fraction of the dark matter velocity 
distribution clS Si function of radius. 

Finally, we note that the unbound fraction of dark matter in halos involves an important 
definitional issue. In the analysis of dark matter simulations to recover dark matter halos, 
many halo finders remove unbound particles (e.g., [1, 44-47]). The usual justification for 
this choice is that unbound particles will quickly leave the system, representing a transient 
population which will not contribute to the long-term mass of the halo. However, some halo 
finders do not remove unbound particles (e.g., [28, 48]), so it is interesting to consider how 
much the observed discrepancies between halo finders (see, e.g., [49]) depend on this choice. 
Moreover, the assumed definition for bound particles traditionally uses the Newtonian escape 
condition (kinetic energy greater than the absolute value of the Newtonian potential energy), 
which has not been appropriately verified in an expanding cosmology. 

These reasons have motivated us to conduct a comprehensive study of the population 
of unbound particles in dark matter halos. In §2, we discuss the dark matter simulation we 
employ as well as the halo finder and conventions used for calculating kinetic and potential 
energies. We present the main results for the population of unbound particles in §3 and 
discuss how these results impact the science considerations above in §4. Finally, we summarize 
our conclusions in §5. Throughout this work, we adopt a flat, ACDM cosmology with main 
parameters = 0.04, Qm = 0.25, Qa = 0.75, and h = 0.7; this is similar to the WMAP7 
best-fit cosmology constraints [50]. 

2 Methods and Conventions 
2.1 Simulation 

We make use of the Consuelo simulation from the Large Suite of Dark Matter Simulations 
(McBride et al, in preparation).^ Consuelo covers a large volume (420 Mpc on a side) 

^LasDamas Project, http://lss.phy.vanderbilt.edu/lasdamas/ 
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with ~2.7 billion particles (1400^), corresponding to a particle mass resolution of 2.7 x 10^ 
Mq . Its large size and resolution make it ideal for particle studies from Milky Way-sized halos 
400 particles) to the largest clusters (IO^^Mq, 400,000 particles). The simulation 
was performed using gadget-2 [51] with collisionless dark matter only, with a force softening 
length of 8 kpc. The assumed simulation conditions were a flat, ACDM cosmology 
(Om = 0.25, Qa = 0.75, h = 0.7, ag = 0.8, and = 1.0), which is similar to the WMAP7 
best-fit cosmology [50]. From this simulation, 100 snapshots were saved from a = 0.075 to 
a = 1.0, spaced at uniform logarithmic intervals in scale factor. 

In a few cases where resolution tests are necessary, we make use of the Esmeralda simu- 
lation, also from the LasDamas Project. This simulation was run with identical parameters 
and software as the Consuelo simulation, with the exception of its box size (640 h-^ Mpc, 
particle count (1250^), mass resolution (1.33 x lO^'' Mq, and force resolution (15 h~^ kpc). 

2.2 Halo Finding and Properties 

Halo finding was performed using the rockstar algorithm, which is a phase-space temporal 
(7D) halo finder designed for high consistency and accuracy of halo properties [1, 52]. The 
ROCKSTAR halo finder locates peaks in the particle phase-space density distribution using a 
locally adaptive phase-space metric. A seed halo is placed at the location of each peak, and 
particles are assigned to the closest seed halo in phase space (see [1] for full details). Particle 
potentials are computed using a tree method (see §2.3), and unbound particles are removed 
before calculating halo properties (spherical over densities, radii, spins, etc.). 

In this analysis, we consider host halos only (i.e., halos whose centers are not within 
the radius of more massive halos). This is because satellite halos (i.e., non-host halos) are 
almost always defined as a collection of self-bound particles within a larger halo: as such, 
they have no unbound particles to analyze. For cases where host / satellite halo relationships 
are ambiguous (such as in major mergers), ROCKSTAR uses the host/satellite relationship at 
the previous snapshot if available. 

We restrict our analysis to halos with more than 400 particles, which have been shown 
to have robustly measurable properties when compared to higher-resolution resimulations 
[53]. Halo virial masses and virial radii are defined according to the spherical overdensity 
criterion of [54]; if a halo includes satellite halos within its virial radius, those are considered 
to contribute towards its virial mass. 

2.3 Kinetic and Potential Energy Calculation 

An individual particle is considered to be bound to a halo if the sum of its physical kinetic en- 
ergy and its potential energy (relative to all other particle positions within the halo radius) is 
negative. Calculating particle kinetic energies is straightforward (|w^ per unit mass), where 
V is the physical velocity of the particle relative to the halo bulk velocity. However, the calcu- 
lation of potential energies in expanding cosmologies poses normalization issues — the usual 
practice of setting the potential to zero at infinity is not well-defined. The standard method 
for bound/unbound classification is to ignore these issues and compute particle potentials in 
a Newtonian metric, so that the potential i;^ of a point mass M (i.e., the Green's function) is 

GM , , 

^{r) = 2.1 

r 

This is the approach taken by, for example, [44-47] and [1]. For completeness, we consider 
both this traditional definition as well as alternate definitions in §3.1, where 6 = for a 
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chosen finite radius Rq; this latter case corresponds to substituting the following Green's 
function: 

GM GM , , 

To calculate particle gravitational potentials, we use a modified Barnes-Hut algorithm [55], 
as detailed in [1]. This method is nearly as accurate as a direct N'^ potential calculation [1]; 
however, as it is a tree-based (©(A^log A^)) algorithm, it is much faster. 

In our analysis, we consider a range of particle binding energies, so it is helpful to define 
the boundedness 6, of a particle as 

\PE\ -KE , , 

(2.3) 



\PE\ 

where KE and PE are the kinetic and potential energies, respectively, of the particle. A 
value of 6 = 1 indicates that the particle has no kinetic energy (is completely bound); a value 
of 6 = indicates that the particle is marginally unbound (kinetic energy equals negative 
potential energy). 

From the virial theorem, one might expect that the average value of b within the virial 
radius of a halo should be 6 = 0.5. However, this is not the case for several reasons: the 
force law in an expanding cosmology is not a simple law (see Eq. (3.2)), not all particle 
orbits stay within the virial radius (see §3.4), and the "virial" radius used here is set via 
a fixed overdensity threshold but cosmological halos are consistently perturbed away from 
virialization by mergers. The combination of these factors makes the average value of b within 
the virial radius of the halos in this study about 0.6 to 0.7. 



3 Results 



We first discuss the statistics of the fraction of unbound particles in halos and explore the 
variation in this fraction with respect to radius, redshift, halo mass, potential definition, halo 
definition and environment in §3.1. Next, we discuss the energy requirements for a particle 
to escape to infinity in an expanding cosmology and compare these to Newtonian energy 
requirements in §3.2. We then consider the particles which actually escape from halos and 
discuss their origins, orbits, and likely future trajectories in §3.3 and §3.4. Finally, we discuss 
the effects of unbound particles on the baryon fraction of clusters in §3.5. 

3.1 Fractional Occurrence of Unbound and Nearly Unbound Particles in Halos 
3.1.1 Average, Radial, and Redshift Dependence 

The definition of an unbound particle is nominally unambiguous: a particle with enough 
kinetic energy to leave its surrounding potential well. However, there can be many different 
choices for what the appropriate potential well should be. Two categories of high-energy 
particles which have particular physical importance are those which can leave the virial 
radius of the halo ("unbound" with (^(i?vir) = 0) and those which can leave the halo and 
never come back ("unbound" with i?i)(oo) = 0, as in Eq. (2.1)). We consider both categories 
of particles here, but we note that most dependencies (on redshift, environment, etc.) are 
very similar for both, so we only show figures for both categories when substantial differences 
occur. In the text and figures, we always specify which potential definition (0(i?vir) = or 
(j){oo) = 0) is being discussed. 
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Figure 1. There are two main definitions of "unbound" — tlie first is tlie ability to escape tlie liost 
halo's virial radius, and the second is the ability to escape to infinity. The top-left panel shows the 
probability distribution of unbound mass fractions according to the first definition (0(i?vir) = 0) as a 
function of halo mass at z = 0; this averages around 30%. The top-right panel shows the probability 
distribution of the unbound fraction which can escape to infinity ((f){oo) = 0) as a function of halo 
mass at z = 0, which averages only about 2-3%. The bottom panel shows the conditional density 
distribution of the latter unbound particle fraction {(j)(oo) = 0) as a function of radius in massive 
< Mvir < IO^'^'^Mq) halos at z = 0. The red dashed line shows the limit of one particle 
per radial bin; effectively, the individual halo sensitivity limit. Because the unbound fraction as a 
whole passes below this limit at small radii, the median unbound fraction (blue dot-dashed line) 
rapidly drops to zero from the floor effect. On the other hand, the average unbound fraction (green 
double dot-dashed line) follows the trend established at higher radii down to a small fraction of 
the virial radius. 




One caveat with this analysis is that, as explained in §3.2, the energy threshold for 
escape to infinity in an expanding cosmology can vary by ±10% at z = (depending on halo 
mass) as compared to the Newtonian escape threshold. To enable better interpretations of 
other results, especially in terms of the "bound" masses reported by halo finders, we keep to 
the traditional Newtonian energy thresholds. In practice, a particle unbound with respect 
to the standard Newtonian escape thresholds for z < 1 will always leave the halo radius and 
will, at minimum, not return for a large number of dynamical times. 
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Figure 2. Top-left panel: The average unbound fraction as a function of mass and radius; smaller halos 
have more unbound particles. For this and other panels, limited resolution prevents calculating unbound 
fractions within 16 h~^ kpc of halo centers. Top-right panel: The unbound fraction as a function of 
radius for alternate halo definitions (i.e., spherical regions enclosing 200 and 500 times the critical density). 
Approximately four times as many particles are unbound if i?5ooc is used as the halo radius instead of Rvir- 
Note that the potential (j> is calculated only using particles within the halo radius; thus, a given particle 
in the halo will appear less bound if a higher overdensity definition is used. Middle-left panel: The 
unbound fraction as a function of radius and redshift for 10^'^ — lO^'^'^Af© halos. The unbound fraction 
appears unchanged to z = 1, but drops off significantly above that redshift. Middle-right panel: The 
unbound fraction for all halos as compared to halos undergoing major mergers. Halos in major mergers have 
surprisingly few extra unbound particles, on average. Bottom-left panel: The average unbound fraction as 
a function of the potential zero-point definition (i.e., the radius beyond which "unbound" particles can travel). 
Most particles have enough kinetic energy to leave the virial radius of the halo, but less than half as many have 
enough kinetic energy to make it to twice the virial radius. Bottom-right panel: The fraction of marginally 
bound particles as a fraction of the boundedness threshold fethresh- For each particle, the boundedness b is 
calculated as its kinetic plus potential energies divided by its potential energy (Eq. (2.3)). A value of fe < 
implies that the particle is unbound; larger values of b (up to b = 1) imply that the particle is more and more 
bound. In this panel, it is clear that the fraction of marginally bound particles is significantly more than the 
fraction of unbound particles (e.g., up to twice as much for particles which are within 10% of being unbound); 
however, it remains extremely difflcult to unbind particles from the central regions of halos. Error bars for all 
panels are calculated using jackknife statistics. 
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Individual host halos can have a wide variety of unbound particle fractions at z = 0, as 
shown in the top panels of figure 1. On average, about 30% of particles in halos at z = 
can escape the halos' virial radii ((/)(i?vir) = 0), whereas only 2-3% can escape to infinity 
{(j){oo) = 0). The width of the distribution in both cases is substantially wider for Milky- 
Way-sized halos (lO^^M© < M^ir < IO^^-^Mq) as compared to massive halos (lO^^M© < 
Mvir < IO^^'^Mq). This is especially pronounced for the <j){oo) = case, where the scatter 
across halos doubles from 0.3 dex in massive halos to 0.6 dex in Milky- Way-sized ones.^ 
We have verified that these probability distributions are identical for our lower-resolution 
simulation. Smaller halos live in a greater diversity of environments; by comparison, massive 
halos will not be found in very isolated environments, nor will they be found near relatively 
much larger halos, which contributes to the larger scatter seen. We address the environmental 
dependence explicitly in §3.1.2. 

We note that some halos in the top panels of Fig. 1 can have very high unbound fractions, 
approaching 100%. These high unbound fractions occur when a smaller halo partially overlaps 
with one or more larger ones. Specifically, all particles within the virial radius of the smaller 
halo will be counted towards its mass, even though those particles may also be within the virial 
radii of the larger halos nearby. These shared particles will have higher velocity dispersion 
from their association with the larger halos, which can lead to very large unbound fractions 
for the smaller halo. 

The bottom panel in figure 1 shows the spread in unbound fractions {(j){oo) = 0) as a 
function of radius for massive halos at z = 0. The significant scatter seen in figure 1 remains 
at all radii. However, at small radii (< 0.1i?vir); the average unbound fraction is so low that 
most halos will have less than a single unbound particle per radial bin. As a consequence, the 
median unbound fraction across halos significantly underestimates the overall trend at small 
radii. For that reason, we express radial unbound fractions in terms of the average across all 
halos, which better follows the trend of unbound fractions established at larger radii. 

The average radial profile of the unbound fraction (i.e., the ratio of unbound to total 
particles in radial bins) is shown in figure 2, in most cases for 4>{oo) = 0. In all cases, the 
unbound fraction increases as a strong power- law with radius; the logarithmic slope ranges 
between 1.5 and 2, steepening towards the outer regions of the halos in question. For a NFW 
[56] halo density profile, this would imply that the physical density of unbound particles 
peaks between one and two times the scale radius. 

The top-left panel in figure 2 shows the mass dependence of the radial dependence of the 
average unbound fraction. As in the top panel of figure 1, smaller halos have larger unbound 
fractions across all radii; at large radii the absolute difference is larger than at small radii. For 
smaller halos in this panel, the radial dependence is truncated at twice the force resolution of 
the simulation (16 h^^ kpc), as the gravitational accuracy of the simulation degrades rapidly 
below that radius. 

The top-right panel in figure 2 shows how the unbound fraction (</>(oo) = 0) changes 
based on the halo mass definition. If the [54] spherical overdensity threshold is replaced by 
another common but higher density definition for massive halos (e.g., 200 or 500 times the 
critical density), then significantly more particles are unbound if the potential is calculated 
only within the halo radius. This is exactly as might be expected: the higher density defi- 
nitions amount to truncating the halo at a smaller radius, reducing the total mass enclosed 
but without reducing the velocity dispersion. The top panel in figure 3 shows for comparison 

^Smaller halos have fewer particles, which introduces a binomial sampling scatter when calculating unbound 
fractions. However, for the smallest halos we consider (IO^^Mq), this scatter is a negligible 0.1 dex. 
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Figure 3. Supplemental panels for figure 2. The top panel shows how the unbound fraction depends 
on the halo mass definition if the potential is defined such that 4>{Rh) = 0, rather than (j){oo) = 0. 
This is equivalent to the fraction of particles which would be expected to have orbits outside the halo 
radius. As compared to the figure for the <^(c») = case in figure 2, the different halo definitions 
result in much more similar unbound fractions. The bottom panels show the evolution of the unbound 
fraction as a function of radius and redshift for two more mass bins (left panel: 10"'^'' — lO^'^'^M©; 
right panel: IO^^Mq — than shown in figure 2. As with the analogous panel in figure 2, 

very little evolution in the unbound fraction is seen to z ~ 1. However, the exact redshift where 
the unbound fraction begins to change does have a slight mass dependence. This may be related to 
the fact that halos which are the dominant sources of mass in their nearby environment have lower 
unbound fractions. 



the fraction of particles which can pass beyond the halo radius {(j){Rh) = 0) for the various 
mass definitions, which is much more stable. 

The middle-left panel in figure 2 shows the redshift dependence of the radial unbound 
fraction, which has apparently no evolution out to z = \. At higher redshifts, however, the 
unbound fraction falls substantially. This trend is preserved across all the mass ranges we 
consider, as shown in the bottom two panels of figure 3. 

The middle-right panel in figure 2 shows how the unbound fraction changes if only halos 
undergoing major mergers are selected. Specifically, a halo is considered to be undergoing 
a major merger if it contains a satellite halo with fmaxjsat ^ O.Guu-iax host) where fmax is 

the 

maximum circular velocity, corresponding approximately to an infalling mass ratio of 0.3:1. 
This panel shows that, for massive halos, the presence of a major merger has very little impact 
on the unbound fraction. However, for less-massive halos (e.g., lO^'^M© < Mvir < IO^'^'^Mq), 
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a major merger can increase the unbound fraction somewhat more (by 1-2% overaU). At 
first glance, it may appear surprising that halos undergoing major mergers do not have 
substantially more unbound particles. However, recalling the definition of the halo mass 
as a spherical overdensity, the mass in an incoming major merger gets added to that of 
the host as it crosses the virial radius. Hence, while the incoming particles may have been 
unbound relative to the potential of the original halo, they are mostly bound relative to the 
combined potential of the incoming and original particles within the virial radius. This effect 
is explicitly shown in §3.1.2. 

The bottom-left panel in figure 2 shows how the unbound fraction would change if the 
zero-point of the Newtonian potential were set at a location other than infinity (see Eq. 
(2.2)). Note that if 4>{Ro) = for a specific choice of Rq, then the unbound particles (6 < 0) 
for this choice of (p are exactly those which have enough kinetic energy that their trajectories 
can reach Rq. For larger values of Ro, the orbital time will increase. As a rough guide, the 
free-fall time of a particle towards a point mass will be v^2i?3(GM)-i, so a particle on a 
radial orbit would spend a time 



outside of the halo. This can be very large even for relatively close distances; particles with 
enough energy to reach LliJ^ir will spend about IGyr outside the halo, and those with enough 
energy to reach 2R^ij- will spend at least lOGyr outside of the halo. Thus, perhaps 20-25% 
of the particles in halos at z = have enough energy to take very long trips outside of the 
virial radius; this is similar to the fraction of particles (25%) which take extended trips out 
of Friends-of-Friends halos [57] . 

Finally, the bottom-right panel in figure 2 shows the radial dependence of marginally 
bound particles. A range of boundedness thresholds (6, see Eq. (2.3)) are considered, from 
b < 0.01 (within 1% or less of being unbound) to b < 0.2 (within 20% or less of being 
unbound). The results in §3.2.2, which account for the effects of halo mass accretion and 
the expansion of the universe, suggest that the appropriate energy threshold for a particle to 
escape to infinity at z = for halos in this mass range (10^^^ to IO^-^-^Mq) is about b ~ 0.05 
at R = and 6 ~ 0.18 at the virial radius. Thus, about 10% of the particles close to the 
virial radius have enough energy to leave the halo and never come back. 

To summarize the most important trends, the unbound fraction is a strong function 
of radius. Regardless of definition, the unbound fraction on galaxy scales (within 0.05i?vir) 
is below 1%, whereas it can reach 10-100% (depending on definition) at the halo radius. 
Major mergers, lower halo masses, and lower redshifts all correlate with increased unbound 
fractions, but these effects are comparatively subdominant. 

3.1.2 Environmental Dependence 

One way to explain the mass dependence of unbound fractions in the top panel of figure 1, 
is through an environmental dependence of the unbound fraction. There are several ways of 
quantifying the local environment; one common method is to use the distance to the nearest 
larger halo. As shown in the top-left panel of figure 4, the unbound fraction for Milky Way- 
sized halos (0(oo) = 0) has a very clear dependence on this measure of environment, with the 
unbound fraction increasing by a factor of 30 over the range of 1-10 Mpc. This dependence 
is even clearer if the distance to the larger halo is scaled by the larger halo's virial radius 
(figure 4, bottom- left panel), strongly suggesting that tidal forces and/or an increase in 




(3.1) 
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Figure 4. Halos in denser environments have larger unbound fractions. Top-left panel: The 
conditional probability distribution of the unbound fraction in 10^^ — lO^^'^A/© halos as a function of 
the distance to the nearest larger halo. Bottom-left panel: Same as top-left panel, except now in 
units of the larger halo's virial radius. This shows that the unbound fraction for smaller halos climbs 
steeply within 3-4i?vir of a larger halo due to the larger background velocity dispersion. Bottom- 
right panel: same as bottom-left panel (the conditional probability distribution of the unbound 
fraction), except for group-scale (10^'^ — 10^^-^ Mq halos. Top-right panel: Unbound fraction in 
10^^ — 1O^^'^M0 host halos as a function of the distance to the nearest major merger or major merger 
candidate, in units of the larger halo's virial radius. For ongoing major mergers, this distance will be 
less than one virial radius; for impending major mergers this distance will be more than one virial 
radius. Once the smaller halo crosses the virial radius of the larger halo, the unbound fraction drops 
swiftly. This is because the mass of the smaller halo is added to that of the larger halo once it crosses 
the virial radius, so that more and more particles are bound to the combined potential of the two 
halos. 



the background velocity dispersion are responsible for the majority of the variation in the 
unbound fraction. Interpreted in the context of dark matter detection experiments, this 
suggests that the high-energy tail of the dark matter velocity distribution function may be 
substantially different for halos which are about to merge into a larger object as compared to 
field halos (i.e., halos which have no nearby neighbors of comparable size). For the Milky Way 
specifically, the nearest larger halo would be associated with M31, which is approximately 
three virial radii away [58]. The average unbound fraction for such halos is 2.8%, practically 
identical to the 2.9% average for all halos in the same mass bin. The Milky Way is thus not 
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Figure 5. The top left panel shows the average unbound fraction for halos as a function of the 
distance to the nearest larger halo, scaled by the larger halo's virial radius (analogous to the bottom 
panels of figure 4). The environmental trend substantially dominates over the mass trend. The top 
right panel shows the probability distribution for the distance to the nearest larger halo in bins of halo 
mass. Unsurprisingly, it is more common for larger halos to be the largest object in their vicinity, and 
therefore not subject to tidal effects which would increase the number of unbound particles. Bottom 
panels: the absolute, as opposed to conditional, probability distribution for the unbound fraction as 
a function of distance to the nearest larger halo (compare to figure 4, bottom panels). 



expected to have a particularly unusual velocity distribution based on its environment. 

In light of the clear environmental dependence of the unbound fraction, it may be 
surprising that halos in major mergers showed comparatively small increases in the unbound 
fraction in figure 2. However, as mentioned previously, particles in one halo which cross the 
virial radius of another will contribute to the virial mass of the second halo. Thus, although 
a major merger is contributing unbound particles, it is also contributing a significant amount 
to the gravitational potential as well. To explore the balance between these two effects, the 
top-right panel of figure 4 shows the unbound fraction in Milky Way-sized halos as a function 
of the distance to the nearest smaller halo with Umax.smaiier > 0.6t>max,iarger- In this panel, 
there is a very clear increase in the unbound fraction as the smaller halo approaches the 
larger one; however, as soon as it begins to contribute a significant amount of mass to the 
larger halo, the unbound fraction rapidly drops back down to nearly the field halo value. At 
the same time, once the smaller halo crosses the virial radius of the larger halo, the majority 
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Figure 6. The evolution of —qH^ with scale factor, relative to Hq. Negative values imply that the 
Hubble flow decelerates escaping particles; escape to infinity is impossible while this is the case (see 
Eq. (3.2)). Positive values of —qH^ imply that the Hubble flow accelerates escaping particles. The 
transition between these states occurs at a = 0.55 for our chosen cosmology. In the far future, —qH^ 
asymptotes to HaHq. 

of the unbound particles which are going to enter the larger halo have already done so. Since 
they are unbound, they will rapidly leave the halo, depleting the high-energy tail of the 
velocity function during the time it takes the smaller halo to sink to the center of the larger 
one. 

It should be noted that the dependence of the average unbound fraction on environment 
for larger halos is very similar to that for Milky Way-size halos (figure 5, top left panel), 
although the field value rises slightly with halo mass. To fully explain the wider scatter 
seen in the probability distribution in figure 1, we note that smaller halos are much more 
likely than larger halos to have a nearby larger neighbor (figure 5, top right panel). The 
environmental impact of a larger neighbor on the unbound fraction is very strong when it 
is within 1-3 virial radii, which leads to a greater scatter overall in Milky Way-sized halos. 
From the bottom panels of figure 4, group-scale halos are more likely than Milky Way-sized 
halos to be the largest halo in their nearby environment, and so are less likely to have high 
unbound fractions, and the same trend holds for cluster-scale halos as well. 

3.2 Energy Thresholds for Particle Ejection 
3.2.1 Analytic Modeling 

As mentioned in §2.3, setting (/'(oo) = in Eq. (2.1) ignores the effects of gravity in an 
expanding universe. Of course, particles in a cosmological halo can easily transfer energy 
and angular momentum through interactions with substructure, so it is perhaps besides the 
point to seek a better analytical solution for particle boundedness. Nonetheless, doing so 
provides a way to better understand how many and which particles could be expected to 
leave the halo and never return in an N-body simulation. 

If one considers a point mass in an expanding cosmology, the effective acceleration at 
experienced by a test particle moving radially with respect to the point mass will be [59]: 

at = (^-^ - qH'ry (3.2) 

where M is the value of the point mass, H is the Hubble expansion rate (da~^), and q is the 
deceleration parameter {—HH~'^ — 1 = —da~^H~'^). Both q and H evolve with time; the 



evolution of the product —qH^ for our chosen cosmology is shown in figure 6. 

At early times (a < 0.55), —qH^ is negative, and it is impossible for the test particle to 
escape to infinity. This is because the deceleration scales linearly with distance: the farther 
the test particle gets, the stronger the deceleration it experiences. However, at later times 
(a > 0.55), —qH^ becomes positive, and it becomes easier and easier to escape. At very late 
times, the expansion of the universe becomes exponential (constant H), implying g — )■ — 1. 
Therefore, the product —qH^ approaches a final value of ^aHq (see also [60]), as may be 
seen from the evolution of in a flat universe: 



(3.3) 



However, the complexity of the time evolution of qH^ means that there is little pedagogical 
value in solving the equation of motion (Eq. (3.2)) exactly. Instead, it is useful to consider 
the late-time limit, i.e., that qH^ is approximately constant. As may be seen in figure 6, this 
is a good approximation even at the present day; qH^ = — 0.63/fo at z = and changes only 
by ~ 10% in the next 3 Gyr. 

Supposing that a test particle has a boundedness 6, the initial velocity will be 



Vq 



'{l-b)2GM 



(3.4) 



and the velocity function of radius will be 



v{r) 



'2GM 



(3.5) 



If there is a solution to the equation v{r > 0) = 0, then the particle will be bound. While it is 
simple to solve the resulting cubic equation in r mathematically, a more instructive solution 
may be obtained intuitively. Within the equivalence radius re, defined as 



GM 



(3.6) 



the gravitational force will dominate and the particle will slow down; however, for r > r^, the 
Hubble acceleration will dominate and the particle will speed up. Hence, if v{re) is positive, 
the particle will escape to infinity; if v{re) is zero or imaginary, then the particle will never 
pass beyond and will therefore be bound. The resulting boundedness condition is therefore 



GM ro 



GM \3 



< 



(3.7) 



This may be readily simplified to obtain 



b > 



2GM 



3^{GMH)i +qH^rl 



(3.8) 



Thus, for any particular choice of initial parameters, it can be quickly verified if the particle 
escapes or not. 
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While Eq. (3.8) only holds for point masses, it is easy to generalize this result. Specifi- 
cally, we may note that the following is a conserved quantity for a test particle: 




(3.9) 



where KE is the physical kinetic energy per unit mass and PE refers to the Newtonian 
potential. Approximating a halo as a spherically-symmetric matter distribution truncated at 
the virial radius (especially valid at late times [61]), we note that Eq. (3.8) can be translated 
using Eq. (3.9) to any radius within the halo as: 



where we have used —GMyiy-R~^l as the Newtonian potential at the halo boundary.^ 

Unlike in a Newtonian potential, the escape speed threshold will depend on the particle's 
trajectory. For non-radial orbits, the particle will have some angular momentum; thus, for 
example, there will be a minimum velocity required at the equivalence radius to conserve 
angular momentum. Specifically, the quantity rvt will be conserved, where vt is the tangential 
component of the particle's velocity. However, obtaining the boundedness condition in this 
case requires solving a quartic equation, so the solution is presented in appendix A in Eq. 



Eqs. (3.10) and (A. 6) are straightforward and simple tests for boundedness. If an 
alternate mass definition is desired, one needs only replace Mvir with the appropriate halo 
mass. If the potential energy is not available for a particle, it may be estimated using Eq. 
(B.3) in appendix B, which gives the Newtonian potential for a particle relative to a [56] halo 
profile. 

It is worthwhile to note that Eq. (3.10) and (A. 6) predict significantly different require- 
ments on boundedness than for the Newtonian case. For example, a test particle at the 
virial radius of a halo at z = could escape on a radial trajectory even if its kinetic energy 
was 34% less in magnitude than its Newtonian potential energy (6 = 0.34)! A particle with 
a purely tangential initial velocity would require slightly over 6% more kinetic energy to 
escape, meaning that any particle at the virial radius with b < 0.29 could escape the halo. 
If the particle were instead at the halo center, the boundedness threshold would depend on 
concentration. Eq. (B.3) would suggest that for a massive halo with concentration c between 
3 — 5, the potential energy at the center would be 3 — 4 times that at the halo radius; this 
would mean that the boundedness threshold would tighten to 6 < 0.084 to 5 < 0.11 to escape 
from the halo center at z = 0. 

Note that, in the late-time limit, the halo mass does not change with time, as all the 
nearby matter has already been accreted [62]. As discussed in §3.2.2, the correction even 
at z = due to the accretion rate of halos is significant, meaning that the true energy 
threshold for escape for particles in a given halo will depend on the individual halo's future 
mass accretion. 

3.2.2 Numerical Simulations 

To investigate the accuracy of the boundedness thresholds derived in the previous section, 
we now turn to numerical simulations of a more realistic model of the halo environment. 

^Note that this assumes that the equivalence radius is beyond the virial radius, which is readily verified. 




(3.10) 



(A.6). 
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Figure 7. The mass accretion history and future for a 10 halo at z = as given by Eqns. 

(3.12)-(3.14). Beyond a scale factor of a = 3, further accretion is practically negligible. 



Specifically, we relax the assumption that the acceleration term —qH^ remains fixed, and 
instead allow it to evolve using Eq. (3.3) and the standard equation for the scale factor a{t) 
in a flat universe: 

a{t) = ^^sinh^ (^^Hoty^^ (3.11) 

We adopt = 0.25, $7a = 0.75, and Hq = 70 km s^^ Mpc^^ for consistency with the 
cosmological simulation we have used. 

As mentioned at the end of §3.2.1, halos are always in the process of accreting matter, 
so the gravitating mass affecting an escaping particle will be changing as well. For these 
numerical simulations, we model halos as having spherical [56] profiles with concentrations 
given by [63] and mass accretion histories accounted for by a slight modification to the 
standard formula from [64]: 

Mvir(a) = Mo exp(z/(Mo, a)) (3.12) 

We allow the extra dependence on redshift on account of the fact that the far future mass 
accretion rates are somewhat less than those predicted by the mass accretion histories in [64] . 
Specifically, based on the mass accretion histories to a = 100 in [62] and the mass accretion 
histories for a < 1 from the merger trees in [52], we adopt the following functional form for 
/(Mo, a): 



/(Mo, a) = 5(a) 



Mo 

-0.122 logio' 



1O6-47M0, 

^ (3.13) 



-0.328 exp (logiOiQii^^ 

r 1 if a < 1 

g(a) = < / 9 (3.14) 

|0.5 + 0.5(^) ifa>l ^ 

An example of the mass accretion history (and future) for a lO^^M© halo at z = is shown 
in figure 7. 

With these assumptions, we can more realistically evaluate the trajectory of a particle 
released in a cosmological environment. An important caveat remains, however, which is the 



assumed definition of "unbound." As mentioned in the previous section, no particle with 
finite energy can ever escape to infinity as long as the deceleration parameter q is positive. 
However, high-energy particles can still escape by waiting until such a time as the deceleration 
parameter becomes negative (for a > 0.55), whereupon it becomes possible to leave the halo. 
We therefore define "unbound" as "can escape to infinity at some time in the future." This 
is well-defined but not trivially so: all halos will eventually stop accreting, and qH^ will 
eventually asymptote to —Q\Hq. As such, there will be a well-defined equivalence radius 
(Eq. (3.6)) which will asymptote to a constant multiple of the virial radius as t — t- oo; particles 
which remain within that radius will stay bound, whereas particles which remain outside that 
radius will escape. 

Another important consideration is the gravitational potential which is used to evaluate 
the motion of particles. Halos are widely thought of as "gravitationally bound objects" [49]. 
Yet, halos do not exist in isolation, and there is a considerable extent of nearby matter which 
contributes to the gravitational potential well (see, e.g., [62]). There are thus two reasonable 
interpretations of "boundedness" : one, that a particle cannot escape from the gravitational 
potential of matter assigned to the halo (e.g., matter within the virial radius); or two, that 
a particle cannot escape from the associated cosmological environment of the halo. We 
consider both definitions, the first corresponding to simulations where we truncate the halo 
density profile at the virial radius ("Virial Potential"), and the second corresponding to a 
truncation of the halo density profile where it matches the critical density of the Universe 
("Full Potential"). 

We show results from the numerical simulations for the mass-dependence of the bound- 
edness threshold at z = in the top panel of figure 8. It is apparent that the mass accretion 
rates of halos represent a non-trivial correction to the result of Eq. (3.8). For smaller halos, 
which accrete less rapidly, the boundedness threshold approaches the value 6 = 0.34 derived 
for a constant halo mass and value for qH^. For larger halos, which accrete more rapidly, 
more kinetic energy is necessary to counteract the future increase in the depth of the gravi- 
tational potential well. Also, as may be expected, more kinetic energy is necessary to escape 
from the full cosmological environment of the halo as opposed to the halo itself. 

As shown in the middle panel of figure 8, the boundedness threshold varies substantially 
with radius. At r = 0, the change in the potential energy barrier due to the expansion of the 
Universe is relatively small, so the boundedness threshold approaches the Newtonian limit 
of 6 = 0. At z = 0, the Universe is accelerating, so particles at larger radii require less 
kinetic energy to escape. However, also at larger radii, the gravitational potential calculated 
from matter within the virial radius (which is used in the definition of b) more severely 
underestimates the full potential of all the matter in the halo's environment. This contributes 
to a feature in the boundedness threshold out to i? ~ 1.5i?vir (at z = 0); if 6 were calculated 
using the full potential instead of the virial potential, then the boundedness threshold for 
the full potential would also rise monotonically with radius. At z = 1, the Universe is 
decelerating; moreover, particles at larger radii feel the effects of deceleration more strongly. 
However, particles at larger radii will also eventually feel more effects from acceleration after 
a = 0.6; hence, the variation of the boundedness threshold with radius is not monotonic for 
either the virial potential or full potential case. 

Finally, the bottom panels of figure 8 show how the boundedness threshold varies with 
scale factor. At early times (a < 0.6), the deceleration of the Universe as well as the rapid 
rate of accretion means that it is very difficult to escape from halos. At very early times, this 
is especially severe; at z = 4, for example, the traditional definition of an unbound particle 
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Figure 8. Boundedness thresholds (&, Eq. (2.3)) for particles relative to the potential of the host 
halo ("Virial Potential") and to the potential of the full cosmological environment ("Full Potential"). 
A boundedness of 6 = corresponds to the traditional definition of boundedness, where the current 
kinetic energy equals the current Newtonian potential energy from the associated halo (i.e., considering 
the potential contribution only from particles within the halo's virial radius). The top panels show 
that the boundedness thresholds vary with mass, in contrast to the result suggested by Eq. (3.8), 
because halo accretion rates vary significantly with mass. The middle panels show how boundedness 
thresholds vary with radius; at r = 0, the traditional definition is very close to being accurate, because 
the effective potential is dominated by gravitation from the halo itself. However, at larger radii, the 
effective potential is increasingly dominated by the expansion of the Universe, to the point that 
particles with very little kinetic energy in comparison to their potential energy (6 — >■ 1) can still 
escape. Finally, the bottom panels show how the threshold varies with scale factor; while the 
Universe was decelerating (a < 0.6), it was much more difficult for particles to escape to infinity. 
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Figure 9. Top left panel: Approximately 2-10% of particles in halos above IO^^Mq at z ^ 1 do not 
appear in the z — descendant. This panel shows the conditional probability density of the escape 
fraction as a function of halo mass a,t z — 1. Top right panel: While a significant fraction of particles 
in 2; = 1 halos do not appear in the z — descendant, the fraction of particles which do not show 
up in any z — halo (as opposed to the main descendant) is only slightly less. This plot shows the 
conditional probability density of the loss fraction as a function of halo mass at z = 1. Bottom left 
panel: for particles within 150 kpc (physical) of the halo center, the escape fraction is substantially 
less, especially for more massive halos. This agrees with the previous section's analysis that even 
marginally unbound particles are rare in the center of massive halos (see figure 2). Bottom right 
panel: for most halos, the average radius (at 2; = 1) of escaping particles between z = 1 and z = 
is tightly clustered at around 70% of the virial radius. This panel shows the conditional probability 
distribution for individual halos at z — 1 for the average radius of particles lost between z = 1 and 
z = 0. 



would underestimate the required kinetic energy for escape by a factor of 4 — 6. This would 
imply that it is not advisable to remove "unbound" particles according to the traditional 
definition at high redshift; however, as shown in figures2 and 3, the error caused by doing so 
is small due to the rapidly dropping fraction of unbound particles. 

3.3 Full Cosmological Simulations 

We now turn to an analysis of which particles actually end up escaping in a full cosmological 
simulation. Since it is impossible to test whether a particle escapes to infinity without running 
the simulation to a very late time, we instead consider which particles can escape the virial 



- 18 - 



5 



Average 




Figure 10. For particles which escape halos between z — 1 and z = 0, this plot shows the conditional 
probability density of the time that they spend in/around halos before escaping (see text), relative 
to the free fall time of the last halo they were in. On average, particles which escaped between z = 1 
and z = spent three free fall times in/around halos before they left. 



radius of their host halos. To do so, we track particles in Consuelo between z ~ 1 and z = 
to see which of them actually remain in the halo with which they were associated at z = 1. 
Specifically, for every progenitor halo at z = 1, we can find the descendant halo at z = 
which receives the largest fraction of the progenitor's particles. The progenitor particles 
which do not end up in the descendant may then be assumed to have escaped sometime 
between z = 1 and z = 0. 

Figure 9 shows the results of this analysis; most halos lose about 2-10% of their particles 
between z ~ 1 and z = 0. This is intermediate between the fraction of particles which can 
escape to infinity (2-3%) and the fraction which can pass beyond the halo radius (30%) 
calculated in §3.1. Intuitively, this may be explained by the fact that halos continually grow 
between z ~ 1 and z = 0, so that some of the particles which would be considered as having 
left the z = 1 halo radius would be counted as being within the z = halo radius. Another 
possibility is that particles could be scooped up by passing halos and end up in a different 
descendant at z = 0; however the top-right panel of figure 9 suggests that this happens for 
only 15-20% of the escaped particles, and that the majority still end up not in any halo at 
z = 0. As shown in figure 9, halos which undergo major mergers between z = and z ~ 1 
(defined as cases where two progenitor halos at z = 1 which have a fmax ratio of at least 
1:0.6 share the same z = descendant) do not have substantially higher loss fractions. 

Considering particles only within the inner 5% of the virial radius (figure 9, bottom 
left panel), the escape fraction is substantially smaller and largely similar to the unbound 
particle fraction within the same radius. In addition, the average radius of particles which 
have escaped between z ~ 1 and z = is roughly 70% of the virial radius; this is significantly 
larger than the average particle radius, which is roughly 40-50% of the virial radius for halos 
in figure 9. 

3.4 Origins of Ejected Particles 

We can gain additional information about escaping particles by considering their origins as 
well as the amount of time they spend in and around halos before being ejected. For this 
analysis, we follow approximately 40 million particles in a subvolume of Consuelo (1/64 of 
the total volume) at every timestep from a = 0.075 to a = 1 and tag particles as either: 
(a) not associated with a halo, (b) smoothly accreted onto its current halo, or (c) initially 
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Figure 11. The origins of most escaped particles are from mergers. Left panel: The fraction of 
escaped particles (between z = 1 and z = 0) which arrived in mergers as compared to the fraction 
which arrived via smooth accretion. Right panel: The fraction of particles within z = and z = 1 
halos which entered the halo via smooth accretion and via mergers. 




Boundedness of Escaping Particles 



Figure 12. The distribution of boundedness levels for particles which leave halos between z = 1 
and z = 0, for the assumption that 0(oo) = 0. Most particles are bound enough that they will never 
escape the halo; their orbit simply takes them outside of the virial radius. However, a small fraction 
of escaping particles have enough energy to never return, even in the far future. 

belonging to a different halo and merging into its current halo. We also tag major mergers 
and cases where the particle leaves the virial radius of its current halo and either returns at 
a later time or enters into a new halo. 

Relevant to the previous section, we limit our analysis to particles which enter halos 
before z = 1 and leave at some point before (but do not return by) z = 0. For each particle, 
we calculate the time it spent in/around halos (i.e., the last time it appeared in a halo minus 
the first time it appeared in a halo) relative to the free fall time of the last halo it was 
in {y^2R^(GM)~^). Figure 10 shows the conditional probability distribution of this time 
for escaped particles in this analysis as a function of the mass of the last halo they were 
in. Regardless of host halo mass, escaped particles spent on average three free fall times 
in/around halos before escaping, corresponding to slightly less than a full orbital period. We 
have verified that this does not depend on whether the last halo the particle was in was 
undergoing a major merger. 
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Figure 13. High-energy incoming baryons will collide with existing gas and remain inside halos, 
whereas high-energy incoming dark matter particles will escape. This will enhance the baryon fraction 
over the cosmic average; this figure shows the relative increase in the baryon fraction as a function of 
halo mass and redshift due to this effect. For small halos, baryonic feedback effects (such as supernovae 
and stellar winds) overwhelm the importance of this effect and decrease the baryon fraction below 
the cosmic average. However, this effect may be more important on the outskirts of massive halos, 
where feedback effects will be less important. Errors shown are jackknife uncertainties in the average 
across halos at a given mass and redshift. 



As shown in the left panel of figure 11, the vast majority of escaping particles came in 
via mergers, as opposed to smooth accretion. The exact fraction is necessarily resolution- 
dependent; we have verified, however, the expected behavior that lowering the resolution of 
the simulation results in a larger smooth accretion fraction. By comparing to the origins 
of particles within the virial radius for halos at z = and z = 1 (figure 11, right panel), 
particles which came in via mergers are much more likely to escape the halo as compared to 
particles which came in via smooth accretion. 

Taken together with the results in previous sections, these results suggest that mergers 
are indeed responsible for ejecting most particles from halos. Although the ejection efficiency 
must decrease for extremely minor mergers, major mergers do not appear to be any more 
effective in ejecting particles than moderate and minor mergers. Finally, we note that the 
majority of escaped particles will eventually return. We show the boundedness distribution 
for escaped particles in figure 12; by comparison the escape thresholds in figure 8, it is clear 
that most of the escaped particles do not have sufficient energy to leave the halo's potential 
well completely. 

3.5 The Effect of Unbound Particles on the Baryon Fraction 

As mentioned in the introduction, high-energy dark matter particles can freely pass through 
dark matter halos, whereas high-energy baryons will shock against the hot halo gas and 
remain inside the halo. This effect will be balanced by galaxy outflows and other baryonic 
feedback effects which will tend to drive baryons out of the halo. These latter effects are 
especially important for lower-mass halos {Mh < IO^^Mq), where the energy transfer in 
feedback effects can be comparable to the binding energy of the baryons. However, in high- 



- 21 - 



mass clusters (M^ > IO^^Mq), the binding energy is large enough that it is very difficult for 
baryons to escape the halo once they are accreted. 

To estimate the effects due to the different collisionalities of baryons and dark matter, 
we evaluate a toy model on all particles in the Consuelo simulation. For dark matter particles 
which are not contained in a halo, we assign the associated baryons (i.e., the cosmic baryon 
fraction times the particle mass) to the last halos in which the particles appeared, if any. 
Then, for all halos, we evaluate the ratio of assigned baryon mass to halo mass to calculate 
the relative increase in the baryon fraction from the collisionality effect. 

This model serves to put an upper limit on the baryon fraction increase from this 
effect; predictions for the relative increase in the baryon fraction as a function of halo mass 
and redshift are shown in figure 13. There is a strong trend towards larger increases in 
the baryon fraction for smaller halos, which may be partially due to environmental effects 
(§3.1.2). As shown in figure 5, smaller halos are less likely to be the largest halo in their 
nearby environment. A nearby, larger neighbor will increase the nearby particle velocity 
dispersion as well as the relative impact velocity of merging satellites — increasing the fraction 
of unbound particles. In addition, because of tidal forces from the larger neighbor, a dark 
matter particle which leaves the immediate vicinity of the halo is unlikely to return. The 
relative change in the baryon fraction also increases at later redshifts; this may be related 
to the smaller fraction of unbound particles at high redshifts (figures 2 and 3) as well as the 
increasing ease of escape at later times due to the acceleration of the universe (figure 8). 

Overall, at z = 0, the relative increase in the baryon fraction is between 15-30%, in 
qualitative agreement with the 25% implied by the results of [57] for Friends-of-Friends 
halos."* For very massive halos, this would suggest an average baryon fraction of nearly 0.19 
within the virial radius at z = 0, instead of the cosmic baryon fraction of 0.16 for our assumed 
cosmology. For smaller halos, we cannot make concrete predictions due to the importance of 
baryonic feedback effects. 

4 Discussion and Implications 

From the perspective of individual particles, the virial radius of a halo is not special; as we 
and others [57] have found, roughly a quarter of the mass of dark matter halos at z = has 
entered and left through this radius. Instead, the most important radius for the individual 
trajectory of particles is the equivalence radius (Eq. (3.6)), beyond which the Hubble flow will 
accelerate the particle away from the gravitational pull of the halo. At early times (a < 0.55), 
the deceleration of the universe makes it impossible for particles to completely escape from 
halos. However, at late times, the acceleration of the universe brings the equivalence radius 
from infinity down to roughly four times the virial radius. This leads to the creation of 
"islands" at late times (a ^> 1), whereby particles (and, indeed, other halos) which are not 
within the equivalence radius are expanded away at rates that gravity cannot overcome (see 
also [60, 62]). 

We have shown that the fraction of unbound particles (regardless of definition) is a 
strong function of distance to the halo center as well as the host halo's environment (§3.1 
and 3.1.2). Perhaps surprisingly, major mergers do not contribute substantially on average 
to the unbound fraction of particles. A major merger at the virial radius causes a substantial 
increase in the unbound fraction, but as the merger passes within the virial radius, its mass 

*While [57] did not make the association to baryon fraction increases, their calculation of total halo masses 
is mathematically identical to that of our toy model. 



- 22 - 



also contributes to the gravitational potential, which raises the kinetic energy necessary for 
a particle to be considered unbound (§3.1.2). The strong environmental dependence that we 
find for the unbound fraction — which is linked to the high-energy tail velocity distribution 
function within halos — is relevant for interpreting dark matter detection experiments. How- 
ever, we have verified that the nearby environment of the Milky Way (i.e., the presence of 
M31) barely changes its expected unbound fraction from the overall average for halos of the 
same mass. 

In terms of precision calibrations of the dark matter halo mass function, we have noted 
that the cosmic baryon fraction will not correspond to the baryon fraction in halos even at 
very large radii (§3.5). This is because the ability of dark matter to freely pass beyond the 
virial radius of a halo is not shared by baryons. As the orbits of many dark matter particles 
take them beyond the virial radius (about a quarter of the mass that has ever been within a 
halo at z = 0) , the baryons left behind will contribute to an excess beyond the cosmic baryon 
fraction in the halo. 

This contribution is mass-dependent, but it is difficult to make concrete predictions for 
the baryon fraction for halos smaller than W^^Mq due to the effects of stellar and quasar 
feedback being important for these halos. However, for high mass halos > W^'^Mq, we 
find a redshift-dependent excess in the baryon fraction on the order of 5-15% (i.e., an observed 
baryon fraction of 0.17-0.19), as discussed in §3.5. This redshift dependence is important 
for observations of clusters, as a changing baryon fraction might otherwise be interpreted as 
an evolving dark energy parameter [65]. Encouragingly, similar trends for how the excess 
baryon fraction evolves with halo mass and redshift are observed in recent hydrodynamical 
simulations [34], although the overall normalization depends on the simulation code used. 

Finally, regarding intergalactic supernovae, we find that it is difficult for central dark 
matter particles to acquire enough energy to leave their hosts (§3.3). Roughly 2-10% of 
particles within halos do so between z = 1 and z = 0, and the fraction drops by an order of 
magnitude for particles within the central 5% of the virial radius. Thus, depending on the 
host halo mass, the fraction of stars (i.e., supernovae progenitors) which receive enough of a 
kick in halo-halo mergers to leave the halo is on the order of a percent or less. Observations 
of hypervelocity supernovae above this fraction should correspond either to stars ejected via 
multi-body stellar interactions or to observational biases. 

5 Conclusions 

We have studied the properties of unbound particles as well as particles which have escaped 
from halos. Our main findings are as follows: 

1. On average, at z = 0, 2-3% of particles in halos can escape to infinity, and 30% have 
enough energy to pass beyond the virial radius. 

2. The fraction of unbound particles in halos is a steep function of distance from the 
halo center as well as environment (specifically, distance to the nearest larger halo). 
Comparatively, it is a weak function of halo mass and redshift. 

3. Based on its halo mass, about ~ 2.9% of the Milky Way's dark matter is expected 
to be able to escape to infinity. The Milky Way itself appears to be typical in this 
sense; if priors about the distance between the Milky Way and Andromeda (M31) are 
included, this fraction changes negligibly. Priors on the distance to larger objects are 
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likely sufficient for determining the unbound fraction of particles contributing to the 
velocity distribution function of dark matter in halos, with relevance for interpretations 
of direct detection measurements. 

4. In simulations, 2-10% of particles within dark matter halos at z = 1 escape and are 
not found within any other halo by z = 0. The vast majority of these escaped particles 
came in from merging satellites, as opposed to smooth accretion. Most of these escaped 
particles do not have enough energy to permanently leave and will return at some time 
in the future within the halo's virial radius. 

5. The baryon fraction in massive (IO^'^-IO-'^^Mq) halos may be boosted within the virial 
radius to 0.19 at z = on account of unbound dark matter particles. This boost is 
redshift-dependent, leading to an expected baryon fraction of 0.175 at z = 1 for our 
assumed cosmology. The standard ACDM cosmological model thus generically predicts 
an evolving baryon fraction with redshift. 
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A The Boundedness of a Test Pcirticle For Arbitrciry Trajectories 

Knowing that the angular momentum rvt will be conserved, we can write down the radial 
velocity Vr of a particle relative to a point mass as a function of its radius r: 




(A.l) 



where vtfl is the original tangential velocity of the particle. For a halo, we can do the same 
outside of the virial radius, where we can write 

vr{r) = - 2Eh - qH'^r^ - (^^^^,0)' (A.2) 

If the particle is unbound, then Vr{r > 0) = can have at most one solution at r < ro; a 
bound particle will have either a solution for r > ro or, in the case of a circular orbit, r = ro 
will be a solution and also a local maximum of v^. One option for distinguishing these cases 
is to search for the zeros of the following quartic equation: 

- gifV^ - 2Ehr^ + 2GMr - r'^vlo = (A.3) 

A somewhat simpler option is to search for the minimum of v'^ and check to see if it is 
negative; the locations of the extrema of f ^ are given by: 



- 2qH'^r^ - 2GMr + 2r^i;|o = 



(A.4) 



Before solving this latter equation, recall that a particle on the threshold of boundedness 
will be one where v1 just barely reaches zero; thus, the location of the minimum of and 
the location of the root will coincide at the turnaround radius. Hence, we can substitute Eq. 
(A. 4) into Eq. (A. 3) to reduce the equation to a quadratic: 

Ehr^ +'^GMr-rlvl^ = Q (A.5) 
This gives a simple boundedness formula for E^'- 

-^GMr + rMn 
Eh < ^ ^2 (A-6) 

We use Ferrari's method to approach the quartic in Eq. (A. 4). The solution is: 



(A.8) 

(A.9) 
(A.IO) 

K = ^ (A.n) 

Of course, there are three more solutions to the quartic; the full solution set is 



y = 






R = 


j + Vj^- 


-K3 


J = 




».6 
' e 
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r=2(±.v/!/±.^-.±13^j (A.12) 

where the ±1 correspond to identical signs, and the ±2 is independent. However, the solution 
in Eq. (A.7) is the only one which corresponds to r = rg (the equivalence radius) for radial 
motion {vtfl = 0). 

These equations have been verified as correct in gravitational simulations of a particle 
escaping a halo. Admittedly, they are quite inelegant, and in some cases, a simpler expression 
may be desired. A useful approximation is to require that angular momentum be conserved 
at the equivalence radius; this results in the following condition for vt^ > 0: 

^,.<-^(«GM™)l + lrKo(ilg)' (A.13) 

For a particle with purely tangential velocity at the virial radius of a halo, this gives the 
correct bound for E^ to within about 1%. 

B The Newtonian Potential of an NFW Halo 

Given an NFW halo [56] with characteristic density p and scale radius r^, the enclosed mass 
at a radius r is 



M(r) = 4TTpr' 



In 



r „ + r \ r 



r + To 



(B.l) 
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Thus, the change in potential from radius ri to radius r2 is 



r2 



GM{r] 



dr 



In 



Ts+r 



(B.2) 



If the halo is considered truncated at the virial radius, the potential outside will be that of a 
sphere with the virial mass. Thus, the combined Newtonian potential as a function of radius 
is: 

r > i?vir 



dGKU 



In 



r 



Rviv 



r<R^ 



where d is the following constant: 



(B.3) 



In 



Rv 



and noting the following limit applies for r — > 0: 



lim 

r->-0 



1 



In 



Vs + r 



(B.4) 



(B.5) 
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